Intrinsic Heating and Cooling in Adiabatic Processes for Bosons in Optical Lattices 
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We show that by raising the lattice "adiabatically" as in all current optical lattice experiments 
on bosons, even though the temperature may decrease initially, it will eventually rise linearly with 
lattice height, taking the system farther away from quantum degeneracy. This increase is intrinsic 
and is caused by the adiabatic compression during the raising of red detuned lattices. However, 
the origin of this heating, which is entirely different from that in the bulk, also shows that one can 
reverse the temperature rise to reach quantum degeneracy by adiabatic expansion, which can be 
achieved by a variety of methods. 
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At present, there are worldwide interests to emulate 
strongly correlated electronic systems using cold atoms 
in optical lattices. Should these efforts be successful, 
we shall have a whole host of new methods for study- 
ing strongly correlated systems which allow one to vary 
density, interaction, and dimensionality with great ease. 
Typically, experiments are performed in the regime where 
the quantum gas is in a tight binding band. A necessary 
condition for reaching the strongly correlated regime is 
that the lattice gas must achieve quantum degeneracy 
in this narrow band. This can be very challenging, as 
the temperature for quantum degeneracy can be several 
orders of magnitude lower than that in the bulk [H. 

In current experiments on bosons, one typically starts 
with a magnetic trap and then turns on a lattice adia- 
batically to bring the system into the Mott regime. How 
temperature changes in this process is a question of key 
interest. While there are theoretical studies 0, there 
are not yet analytic understanding of the cooling power 
of this process. For a non-interacting gas in a tight bind- 
ing band, the only energy scale is the hopping integral t, 
which decreases exponentially with lattice height V a . By 
dimensional analysis, the entropy density of a quantum 
gas in an infinite lattice must be a function of T/tQ, 
where T is the temperature. At first sight, this seems 
to furnish a powerful cooling scheme, for t will drop ex- 
ponentially with increasing V , which will force T to do 
the same to maintain adiabaticity. This, however, does 
not work in practice. The entropy density is not only a 
function of T/t but also U/t, where U is the on-site re- 
pulsion energy. As both t and T decrease, U/t becomes 
so important that it overwhelms the T/t contribution. 
Note that interaction effect is not only important in the 
Mott regime, but also in the superfluid phases close to 
the quantum critical point where condensate depletion 
becomes severe. Thus, even on the superfluid side, the 
ratio T/t can not be constant for sufficiently large V a . 

As a matter of fact, in an infinite lattice, the temper- 
ature of Bose gas must rise when it is brought deeper 
into the Mott regime. This is because the energy gap in 
the Mott phase is given by U, which increases with lat- 
tice height Vo^jfy. Raising the lattice will therefore make 
it harder to generate excitations. The only way to keep 



the entropy constant is then to raise the temperature to 
counter the rising excitation energy, which unfortunately 
drives the system away from quantum degeneracy. We 
shall call this "temperature runaway" . 

The presence of an harmonic trap, however, has pro- 
found effect on the adiabatic processes, to the point that 
the bulk effects mentioned above are completely irrele- 
vant. In a trap, a lattice Bose gas will have a "wedding 
cake" structure consisting of concentric regions of Mott 
phases. Between neighboring Mott regions is a shell of 
mobile atoms (referred to as a "conducting" shell) . Since 
boson number fluctuates in these shells, they have much 
higher entropy density than the Mott regions at temper- 
atures T < U, and are the sources of entropy at low 
temperatures. As we shall see, the adiabatic processes in 
current experiments produces an adiabatic compression 
on these conducting shells which again cause a temper- 
ature runaway. However, by understanding how entropy 
distributes in the system, it is possible to reverse this 
"runaway" by a simple process, i.e. adiabatic expansion. 
In this way, one can reach down to the degenerate tem- 
perature of the order of t even for deep lattices. 

(A). The relevant energy scales: The po- 
tential of an infinite optical lattice is V L (x) = 
V Q Xa=i 2 3 sm 2 (Trxi/d) 7 where d is the lattice spacing. In 
deep lattices, bosons will reside in the lowest band and 
the system is described by the boson Hubbard model. 
In grand canonical form, it is K, = H t + Hjj — fiN 
= K, + H t , where H t = -iE(a,H') = 
J2r Un R (n R — l)/2, ]C = TLu — uN , a R creates a boson 
at site R, n R = a R a R , and fi is the chemical potential. 
Both the hopping integral t and on-site energy U can 
be calculated from lattice height V and recoil energy 
E R = {Trh) 2 /{2Md 2 )\^. To give an idea of the energy 
scales in the Mott regime, we show the values Eg, U, 
t and t 2 /U in Table I for 87 Rb (with scattering length 
a s = 5A5nm) in a lattice with d = 425nm for various 
lattice height, where Eg is the lowest band gap in the 
density of state, and t 2 /U is the scale for virtual hop- 
ping. 

We see from the Table I that for V /Er > 10, we have 
U >> t >> t 2 /U. Since the superfluid-insulator transi- 



2 



tion is estimated to be V /E R = 12-13 for 87 Rb0, the 
condition U » t » t 2 /U \s well satisfied in the region 
of transition, and is strongly enforced at higher lattice 
height. Table 1 shows the great challenge of reaching 
quantum degeneracy (i.e. the scale t^) even for V as 
low as 15-Er, and the even greater challenge of reaching 
temperatures — t 2 /U. 
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(B) . Adiabatic compression: In current experi- 
ments, the optical lattices are constructed from red de- 
tuned lasers where atoms are sitting in the region where 
the laser density is high. The lattice is switched on adi- 
abatically in a magnetic trap V m {r) = Muj 2 r 2 j2, with 
frequency lj . However, due to the Gaussian profile of 
the laser beam, the laser itself will also produce a confin- 
ing harmonic potential. When these two potentials are 
properly aligned, the frequency w of the total harmonic 
trap V(r) = Muj 2 r 2 /2 is@ 

lo 2 =u 2 + W /{Mw 2 ), (1) 

where w is the waist of the laser beam. Thus, as the 
lattice is switched on, it also provides an adiabatic com- 
pression. In fact, for strong lattices, the laser contribu- 
tion can dominate over that of the magnetic trap. For 
example, for V = 12Er, Er = h X 3.2kHz, w = 130/^to, 
and loo — 2tt x 15Hz, the second term in eq.([T]) is already 
10 times larger than the first. As we shall see, this adia- 
batic compression is the source of temperature runaway 
in a trap. 

(C) . Density profile and entropy distribution: 

We shall focus on the case U » t, and temperatures 
U >> T > 0. Let us distinguish the cases T > t, t > T, 
and T = 0. (T has the dimension of energy). 

(CI) T = 0: If t were zero, all sites are decoupled. The 
eigenstates on each site are number states. The ground 
state then has = m, (m = 0, 1,2..), when a lies in 

the interval (to — 1)U < \i < mU. In a trap, the density 
profile can be obtained using local density approximation 
(LDA) by replacing /i with /z(r) = fj, — V(r). This leads 
to a "wedding cake" structure, (see figure 1), with sharp 
steps located at 

R m = ^2{ f i-mU)/{Mu J 2 ), (2) 

and u determined by the number constraint, 

N= /drn(r) = ^£0R m /d) 3 (3) 

J 6 m 

The m-sum are restricted to < m < fJi/U. 
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FIG. 1: Density distribution in a harmonic trap. The density 
profile of the steps at T > t is described by the Fermi function 
in eq.@. The system considered has 3.5 x 10 5 87 Rb bosons, 
in a magnetic trap with uj — 2n{15Hz), V /Er = 25, d = 
425nm. 

Since the energy exciting from m to (m + 1) boson 
state is £ m +\ — £m = fnU — u, the number states m and 
m+1 are degenerate at u = mil . This degeneracy will be 
lifted when t ^ 0, which introduces fluctuations between 
neighboring number states. As a result, the sharp steps 
in fig.l is rounded off over a shell of width {AR m )(°) — 
(to + \)t/{Muj 2 R m ). (See (C6) for more discussions). 

It is also important to note that in addition to the 
adiabatic compression, the increase of V will makes the 
Bose gas more repulsive, since U — {V/ Er)°' S8 ~ 140]. 
As a result, /i also increases almost linearly as V (see 
eq.©) and the radii R m only has a weak V dependence. 
In other words, the trap tightening effect and increase 
of repulsion, both caused by the rising V , almost cancel 
each other. This effect can be verified numerically, and 
is important for our later discussions. 

(C2) U >> T » t: In this case, the thermodynamics 
is again dominated by on-site interaction /C Q with H t as 
a perturbation. It is straightforward to show that for 
fi — mU, the number density n{T, a) = (a'a)/d 3 is, 

n(T, M ) = rf- 3 [TO + /] + (..), (4) 

where / stands for f{mU — fi), f[x) = {e x / T + l) -1 , and 
(..) denotes terms of order of ((t/T) 2 , e~ u ' T ) and higher. 

The fact that the density profile of the step is a Fermi 
function with width T means that this width can be used 
as a temperature scale of the system. It is also simple to 
show that the entropy density is 

s{T, fi) = d- 3 [-(1 - /)ln(l - /) - /In/] + (..). (5) 

Eq.([J| shows that the location of the conducting shell 
(at / = 1/2) is still given by (eq.©). Applying the 
Sommerfeld expansion on eq.© shows that the relation 
between n and N is still given by eq.([3|) with corrections 
0(T/U) 2 . The width of the conducting layer AR m is 
now given by A^t(r) ~ T, or 

AR m ~ T/{Mio 2 R m ). (6) 
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FIG. 2: Entropy density as a function of position with the 
same parameters as in fig.l. The solid and dotted curves are 
for T — 5nK and T = InK. The decrease of the peak value 
of s at smaller r is due to the (t/T) 2 and Sommerfeld expan- 
sion correction in eq.©. The inset is the entropy density at 
temperatures T < t discussed in section (C6). 



In figure 2, we have plotted the entropy density as a 
function of position. One sees that the entropy density 
s(r) is concentrated in the conducting layers and is essen- 
tially zero (~ 0(e~ c/ / T )) in the Mott region. The maxi- 
mum value of entropy per site s(r)d 3 is ln2 + 0((t/T) 2 ). 
It occurs at R mi reflecting the degeneracy of m and m+1 
number states at fi = mil. The total entropy of the wed- 
ding cake structure S ca ke = 4-7T J drr 2 s(T, /z(r)) is then ~ 
4 ^E„X,Ai? m ln2/d 3 , or S cake ~ J2 m TR m / (mu 2 d 3 ). 
Using Sommerfeld expansion, and with the (t/T) 2 cor- 
rection, we have [a 
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T 



3 Moj 2 (P 



R m (l + 0((t/T) 2 ,(T/U) 2 )). 

(7) 

(C3) Temperature runaway: Eq. (J7J) is the origin 
of temperature runaway. Since to 2 increases with lat- 
tice height (see eq. (JTJ) ) , T must increase to keep the en- 
tropy constant. The simplest case is a single Mott re- 
gion, which occurs when U/2 > fi > 0. In this case, 
Ra = y/2n/(Mu) 2 ), and N = 4tt R%/ (3d 3 ). The size of 
Mott regime (Rq) is independent of T and to; and eq.© 

becomes S cake (T) = (^) 4 . In the limit of 

large V a , eq.lJTJ) implies S ca k e oc T/V a . In general, there 
are several Mott steps. However, as we have discussed 
at the end of (CI), the sum ^2, m Rm in eq.© depends 
weakly on V . So we again have S ca ke oc T/V , hence the 
same temperature runaway. 

(C4) Reversing the temperature runaway: Since 
the temperature increase is caused by adiabatic compres- 
sion, the simplest way to reverse it is to reduce uj 2 , i.e. 
an adiabatic expansion. Some possible ways are: (a) 
adding a blue detuned laser to generate a repulsive trap 



(with negative curvature 
that eq. (JTJ) becomes to 2 — i 



') of variable strength so 
u\ +8V /mw 2 . (b) Start 



sion can then be achieved by increase V . (c) Turn on a 
(magnetic) anti-trap which make the sign of u> 2 negative. 
This can be achieved by flipping the spin of the boson 
from low field seeking state to high field seeking. This 
process is adiabatic as the flipping of spin does not affect 
the number distribution of the system, (d) Use a laser 
beam with large width w and reduce uj q . The reduction 
of oj 2 will enable one to reach the lowest possible temper- 
atures within the regime U > T > t. As T — ► t, eq.0 
is no longer valid, and the entropy has to be calculated 
differently. (See (C6)). 

(C5) Condition for cooling: Before discussing the 
case T < t, let us illustrate how the initial temperature 
Tjj of a Bose gas in a magnetic trap without lattice is 
related to the temperature Tf of the final wedding cake 
structure. The thermodynamics of the a Bose gas in 
harmonic trap (no lattice) has been studied by Giogini, 
Pitaevskii, and Stringari[9( using Popov approximation. 
For initial temperature Tj < \x , where fi a is the chemical 
potential in the center of the trap, their results imply the 
entropy is 
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with a very tight harmonic trap and use a blue (instead of 



where a = yhj (Mco ) and a s is the s-wave scattering 
length, and A — 10.6. 

The temperatures T{ and Tf are related as Si(Ti, N) = 
S cake (Tf,N; V ), where S cake (T, N; V Q ) is obtained by in- 
verting the number relation N = J n(T,/i(r)) = N(T,/j,) 
for jj, and substituting it into the relation S ca ke — 
J s(T,fi(r)). In figure 3, we have represented Si by a 
black solid line for a system of 3.5 x 10 5 87 Rb bosons, 
lu = 2ir(15Hz), a s — 5.45nm. Tj is chosen to be 13nK. 
The entropies S ca ke of the final states obtained by raising 
a red laser lattice to V /Er = 15 and 30 in the same mag- 
netic trap are represented by a dotted and dotted-dash 
curve respectively. At these lattice heights, the system 
has three Mott layers. The reason that S ca ke deviates 
from the linear dependence on T when T increases is 
because of the corrections to Summerfeld expansion in 
eq.Q and (JSJ). The dashed line with large slope is the 
entropy of a final state with low (total) trap frequency 
(lu = 27r(15)Hz), produced by the expansion schemes 
mentioned in Section (C3). In this case, the system has 
only one Mott layer. 

In the entropy-temperature plot in fig. 3, Ti and Tf 
are connected by a horizontal line. Heating (cooling) 
occurs if S ca ke lies on the right (left) hand side of Tj. 
Temperature runaway is reflected in the fact that the 
slope of S ca ke(T, N;V ) decreases with increasing V , and 
will eventually lie on the right hand side of any T. For 
given T, the critical V above which intrinsic heating 
occurs is given by S^T^N) = S , ca fe e (T 1 l , N; V*). For our 
example, this critical V* is about 25-Er. From fig. 3, it 
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FIG. 3: Illustration of intrinsic heating and cooling: Solid 
curve is Si. The dotted and dotted-dashed line are for 
Vo/Rr = 15 and 30 resp. The system has the same parame- 
ters as that in fig.l. Point (a) denotes the initial temperature 
Ti. The final temperature for V /Er = 15 and 30 are de- 
noted by point (b) and (c). The dashed curve with steepest 
slope is the total entropy for the adiabatic expansion process 
with uj = 2ir(15H z). The final temperature is point (d). For 
Ti > VbnK, considerable heating will occur when raising a 
red laser lattice beyond V = 15Er. Adiabatic expansion, 
however, can cool the system considerably. 



is clear that if the initial temperature is above 20nK , 
raising a red laser lattice will lead to a significant increase 
in temperature. (See also footnote [§1). In contrast, if we 
use a blue laser for the lattice, one can gain an order of 
magnitude or more reduction in temperature. 

(C6) T < t: In this case, hopping is not a pertur- 
bation. We shall focus on the regime /i ~ mU where 
entropy density is highest. This regime is intricate for 
the following reasons. At \i = mU, the ground state of 
the Bose Hubbard model is a superfluid, with transition 
temperature T c is of the order (m + l)t. As the difference 
|/i — mU\ increases, number fluctuations are suppressed 
and T c drops. At some point, the falling T c will reach T 
and drop below it. To calculate the entropy as a func- 
tion of T and /x, one must face the complexity of critical 
phenomenon at some point. While it is possible to cal- 
culate the entropy function using quantum Monte Carlo 
methods, it is desirable to have analytic understanding 
so as to perform quick estimates, which we do below. 

Near /i ~ mU , dominant number fluctuations are be- 
tween the number states m + 1 and m on each site. 
Denoting these two states as a "pseudo-spin" |l/2) and 
| — 1/2), we have or = S^, ol a-p = + m + 1/2; 
and the Hamiltonian K written as[lfj, K. = —h Y^r ^r ~~ 
E<r,r^Sr-Sr'+^J^,] where = (S£,S£), h 
= fi-mU + 0(t 2 /U), J ± = (m+l)t, J z = 0(t 2 /U). Su- 
perfluid order corresponds to (S^) ^ 0. The Mott phase 
corresponds to (Sr.) = ±1/2. 

At h — 0, Jj_ ^ 0, the system at T = is a fer- 
romagnet in the xy-p\a,ne (superfluid). The transition 
temperature is T c ~ J±. As \h\ increases, the order- 
ing spin will tilt away from the xy-plane and develop 



an S z component. At the same time, the magnitude 
of (Sr) is reduced, and so as T c . Using Primakoff- 
Holstein transformation, it is straightforward to work 
out the spin wave spectrum and hence the entropy den- 
sity at temperatures T < t near fi ~ mU, which 



is s, 
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and 



{m+l)t J \ %t(m+l) / 

where C = V3tt 2 /1620. The fact that s increases as 
h moves away from reflects the weakening of superfluid 
order when h ^= 0. 

For \h\ > J±, (or \n — mU\ ^> (m + l)t), the system 
enters the Mott phase. Since T < t ~ J±, the entropy 
density will reduce to the Mott value (which is essen- 
tially zero) over a chemical potential range A(fi—V(r)) ~ 
(m+l)t, or spatial distance AR m ~ (m + l)t/ (Miv 2 R m ). 
(See the inset of fig. 2). The case \h\ ~ J is diffi- 
cult one. Accurate answers will require quantum Monte 
Carlo treatments. However, with the features of en- 
tropy density mentioned above, one can estimate the 
total entropy to be S cafce ~ i'7rJ2 m R 2 n (AR m )s m , or 

Scake ~ 4 Hf MJW Em^TW- ThlS faSter dr0 P ° f 

Scafce with temperature will reduce the cooling power of 
adiabatic expansion. However, this will take place at 
temperature scale so low that it is not visible on the scale 
of fig. 3. 
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